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Abstract 

One of the remaining open issues in the context of the analysis of Active Galactic Nu- 
clei (AGN) is the evidence that nuclear gravitational accretion is often accompanied by a 
concurrent starburst (SB) activity. What is, in this picture, the role played by the obscur- 
ing dust around the nucleus and what do the state of the art AGN torus models have to 
say? Can the IR data provided by Spitzer and Herschel help us in extensively investigate 
both phenomena and, if so, how and with what limitations? In this paper we present our 
contribution to the efforts of answering these questions. We show some of the main re- 
sults coming from a comparative study of various AGN SED modeling approaches, focusing 
mostly on the much-debated issue about the morphology of the dust distribution in the 
toroidal structure surrounding the AGN. We found that the properties of dust in AGN as 
measured by matching observations (be it broad band IR photometry or IR spectra) with 
models, strongly depend on the choice of the dust distribution. Then, we present the spec- 
tral energy distribution (SED) fitting procedure we developed, making make the best use 
of Spitzer and Herschel SPIRE mid- and far-IR observations, to dig into the role played by 
the possible presence of an AGN on the host galaxy's properties. 

1 Introduction 

The central black hole of Active Galactic Nuclei (AGN) is belived to be surrounded by a 
dusty structure, commonly referred to us as "AGN dusty torus", which, according to the 
unified scheme for AGN [Ij plays a crucial role in explaining the difference between type 1 
and type 2 AGN as an orientation effect |24j . The intersection of the line of sight with the 
obscuring dust can give rise to the observed lack of broad lines in the optical spectra, typical 
of type 2 objects, and to a number of differences in the observed SED of AGN at almost all 
wavelengths. 



2 



Modelling IR SED of AGN 



The various radiative transfer models of AGN tori that have been developed in the 
last decades can be divided mainly into two categories: those in which the dust is uniformly 
distributed in an axisymmetric geometry around the central supermassive black hole, cus- 
tomarily referred to as "smooth" models" [22l[Hl[MlEllIZlEIlEnilIj,the first that has been 
developed because more simple to compute and those in which the dust is distributed in 
clumps or clouds surrounding the central source, and hence referred to as "clumpy models" 
\18\ m \T3\ I27j . In the first part of this paper we report the main results coming from a 
systematic comparison (see Feltre et al. 2012 for more details [6]) between two of the most 
widely used models in the literature, each one representative of one of the two classes. 

In order to get a more complete view of the phenomena we also perform a multiwave- 
length study developing an ad hoc SED fitting procedure aimed at reproducing the SEDs of 
the AGN from the optical/UV to the far-IR. Thanks to its multi-component approach this 
tool turned out to be very powerful for the investigation of the interplay between the AGN 
and the host galaxy. The methodology and the first preliminary results are presented in the 
second part of the paper. 



2 AGN Torus models comparison 

In order to investigate the possibility to use the SED fitting technique as a diagnostic for the 
dust distribution we carried out a systematic comparison between two of the most popular 
models in the literature, namely: smooth models of Fritz et al. 2006 |7j, hereafter FOG, 
and clumpy models of Nenkova et al. 2008 |T8| [T9l I20j . hereafter N08. We first conducted 



an automatic comparison (Sec. ^2.2) investigating the main properties of the model SEDs 



and, in a second step, we compared the models with observations considering mid-IR spectra 



information and photometric data (Sec. ^2.3). 



2.1 Grid selection 

The comparison was carried out considering two grids of the aforementioned models, FOG and 
NOB, selected trying to match the values of the model parameters, such as the outer-inner 
radius ratio, the optical depth along the line of sight and the torus aperture. Furthermore, in 
the model-to- model comparison we only considered two extreme inclinations, corresponding 
to type 1 and type 2 objects. Fig. [T] shows region occupied by the two model grids for type 
1 objects and type 2 views. 



2.2 Model-to-model comparison 



In Feltre et al. 2012 [Gj we investigated different properties of the IR SED of the models, that 
is the width and the peak wavelength of the IR SED, the strength of the silicate features at 
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Figure 1: The range of SEDs covered by the smooth (left) and clumpy (right) dust con- 
figurations in the restricted parameter grids. The coverage is shown for both type 1 (filled 
regions) and type 2 (dashed) inclinations 

9.7 and 18 /xm[^ the spectral index in the mid-IR a/i^and the luminosity at 12.3 /im. Our 
findings can be summarised as follows: 

- S9.7 and Sis of both F06 and N08 models are shown in the feature-feature diagram in Fig. 
[2] for both type 1 (left) and type 2 (right) views. In particular, for type 1 views F06 and 
NOB models occupy different region of the diagram. This difference arises from the fact that 
the two models use different silicate absorption coefficients: Laor & Draine 1993 |15) and 
Ossenkopf 1992 [21] for F06 and N08, respectively; 

- the distribution of the width of the IR SEDs, Wm, is narrower for N08 models, for both type 
1 and 2 views, despite the presence of a range of values common with the Wm distribution 
of FOG models; 

- the distribution of am, especially for type 1 views, is very different between F06 and N08 
models and the common range of values is very narrow. This difference is mostly likely due 
to the two different power laws used to model the primary source and lack of a very hot 
component in the clumpy models (see e.g. [3j). 



^Defined as the logarithm of the ratio between the flux F measured within the Une profile over the continuum 
flux Fc at such wavelength, i.e. 



where \m is the wavelength at which the feature's strength is an extremum with a value in the interval between 
8.5 and 11.5 /xm for Sg.? and between 17.0 and 19.5 /.im for Sig-lfor the computation of Fc see Sirocky et al. 
2008 128]. FOG and N08 models consider different absorption coefficients that peak at different wavelengths, 
but for simplicity we call the strength of the silicate feature around 9.7 pm 5*9.7, irrespective of A^. 



5' = ln(F(A™)/Fe(A„)) 



(1) 



■^Defined as 



OlIR = 



logio(F4.5) -logio(F3.6) 
loglo(^4.5) - logio(A3.6) ' 



(2) 
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Figure 2: 5i8 versus Sg.y for smooth (green squares) and clumpy (red triangles) models, and 
89,7 and Sis distributions for type 1 (left) and type 2 (right) views. [6j 




Figure 3: Wjr and am, for type 1 (left column) and type 2 (right column) views, for smooth 
(continuous green lines) and clumpy (dashed red lines) models. [6] 
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Figure 4: Object SDSS_13635.0 (RA = 10:34:13.89(hms), Dec = +58:52:52.8(dms)) is repro- 
duced by a single component, namely that of a type 1 AGN FOG model (left) and NOB model 
(right). [6] 
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Figure 5: Example of the SED fitting of the IRS spectrum of AGN24 [SH] for F06 models 
(top panel - dark green line) and N08 models (bottom panel - red line).|6j 

2.3 Comparison with observations 

The second step is to compare the models with observations. We firstly consider a sample of 
278 spectroscopically confirmed SDSS type 1 quasars with redshift spanning 0.06 < z < 5.2 
with photometry coming from SDSS, 2MASS (whenever available) and SWIRE (Spitzer Wide- 
area InfraRed Extragalactic Survey |16j ) . This sample has been presented in Hatziminaoglou 
et al. 2008 [9]. We also consider a sample of 160 objects with mid-IR spectra which was a 
part of a collection of IRS spectra available in the literature (covering the wavelength range 
from ~ 5 to ~ 40 /im) presented in Hernan-Caballero & Hatziminaoglou 2011 jl2j. As can be 
seen in the examples reported in Figs. [4] and [5] both "smooth" and "clumpy" models provide 
equally good fits when photometric data or IRS spectra are considered. 



3 Multiwavelength emission of AGN 

The availability of a large quantity of photometric data points available requires a customised 
and versatile tool to investigate the multi-wavelength properties of extragalactic sources. We 
present here a multi-component multi-wavelength SED fitting procedure accounting for the 
different emission mechanisms acting from the UV to the far-IR. 

3.1 SED fitting procedure 

We have very recently made further improvements to a fully automatic fitting procedure, 
described for the first time in Hatziminaoglou et al. (2008) [9j and aimed to reproduce data 
from the optical/UV to the far-IR/submm. The procedure has been already applied in various 
published works [101 El Ell [IS [23] and extensively tested with all the different options in 
the context of the ifersc/ie/ Multi-tiered Extragalactic Survey (HerMES; PI S. Oliver). Three 
emission components (see also Fig. pi) are considered: 
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Figure 6: Example of a best fit to a low (top) and a high (bottom) redshift object. The 
IRS spectrum (light blue; see also the insert figure) and photometric data (red symbols) 
are reproduced using the three emission components: SSPs (dotted dark green), AGN torus 
(dashed dark blue line) and starburst (dashed light green line), all together giving the total 
model emission (in black). 



- simple stellar population (SSPs) models, built using the Padova evolutionary tracks [2], 
account for the optical-UV emission; 

- AGN torus models, either smooth (FOG) or clumpy (N08), to account for most of the mid-IR 
emission, due to hot dust heated by the presence of an AGN; 

- empirical starburst galaxy templates to reproduce the bulk of the far-IR emission dure to 
the presence of cold dust. 

A recently introduced improvement is the possibility to fit simultaneously mid-IR 
Spitzer/IRS spectra and broad band photometry. IRS spectra provide constrains on the 
torus models, especially when a silicate feature is present in emission or absorption, while 
Herschel photometry, especially in the long (SPIRE) wavelengths, allows us to constrain the 
peak and width of the cold dust emission. This can be seen in the illustrative examples of fit 
reported in Fig. [6| Furthermore, in order to compute the properties of the cold dust, such 
as the temperature and the mass we reproduced, in a second step, the far-IR emission with 
a grid of modified block body (with fixe value of /3 = 2). This is a very powerful tool as it 
allows us to derive various physical properties for each of the components, such as the AGN 
accretion luminosity, the AGN and the starburst contributions to the total IR luminosity, the 
optical depth of the hot dust, the size and the mass of the torus, the stellar mass of the host 
galaxy and the star formation rate. 



3.2 Study of the AGN impact on a large sample 

We are now applying the full procedure to a large sample of ~500 galaxies with multi- 
wavelngth photometry (from SDSS to Herschel) and IRS spectra available. The aim of this 
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work is to study the impact of a possible presence of an AGN on the host galaxy's properties, 
analyzing the relative contribution of AGN to the total emission, the star formation rate in 
connection to the AGN and the mass of the hot (AGN) and cold (starburst) dust components. 

Preliminary results shows a consistency, despite a non-negligible scatter, between dif- 
ferent estimates of the star formation rate (SFR) (e.g. from the IR luminosity, from PAH 
and from the fit to the optical data) and their independence from the AGN contribution to 
the total far-IR luminosity. Furthermore, no correlation has been found between the masses 
of the hot and cold dust and between the AGN accretion luminosity and the temperature of 
the cold dust (coming from the fit with the modified black body). 

4 Discussion and Conclusions 

The conclusions we obtained from our results can be summarized as follows: 

- the differences in the properties between FOG and N08 models arise mostly from the models 
assumptions (e.g. primary source, dust chemical composition) and not from the dust mor- 
phology with the consequence that the hot (AGN) dust properties, such as the size and the 
mass of the torus, strongly depend on the choice of the torus model; 

- multi-wavelength SED fitting alone does not allow to distinguish between the two dust 
distributions; 

- the independence of the computed SFRs on the AGN contribution to the total far-IR 
luminosity indicates that the AGN component does not infiuence the star formation process; 

- the absence of correlations found between the properties of the hot (AGN) and cold (SB) 
dust is an indication that the two dust heating mechanisms act on different scales. Indeed 
the central engine heats the dusty torus within few tens of pc at most, while the young stars 
in star forming regions, which are mainly responsible for cold dust heating, extends out to 
kpc scales. 
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